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Abstract 

Traditional cutoff regularization scliemes of the Nambu-Jona-Lasinio model limit the applica- 
bility of the model to energy-momentum scales much below the value of the regularizing cutoff. In 
particular, the model cannot be used to study quark matter with Fermi momenta larger than the 
cutoff. In the present work an extension of the model to high temperatures and densities recently 
proposed by Casalbuoni, Gatto, Nardulli, and Ruggieri is used in connection with an implicit regu- 
larization scheme. This is done by making use of scaling relations of the divergent one-loop integrals 
that relate these integrals at different energy- momentum scales. Fixing the pion decay constant at 
the chiral symmetry breaking scale in the vacuum, the scaling relations predict a running coupling 
constant that decreases as the regularization scale increases, implementing in a schematic way the 
property of asymptotic freedom of quantum chromo dynamics. If the regularization scale is allowed 
to increase with density and temperature, the coupling will decrease with density and temperature, 
extending in this way the applicability of the model to high densities and temperatures. These 
results are obtained without specifying an explicit regularization. As an illustration of the formal- 
ism, numerical results are obtained for the finite density and finite temperature quark condensate, 
and to the problem of color superconductivity at high quark densities and finite temperature. 

PACS numbers: 12.39.-x, 24.85+p, 21.65. +f 



I. INTRODUCTION 

The Nambu-Jona-Lasinio (NJL) [l| model has been the prototype model for studying 
chiral symmetry restoration in hadronic matter at finite baryon densities pn and finite tem- 
peratures T. Since the earlier applications of the model at high ps and T [2, l3| , an extensive 



and important body of work has 



)een done in this direction - for reviews and a comprehen- 



sive list of references see Refs. [4, |5|, la, [Zl, Isl, la, [lO|]. Because the model is nonrenormalizable, 
the high momentum part of the model has to be regularized in a phenomenological way. 
The common practice has been to regularize the divergent loop amplitudes with a three- 
dimensional momentum cutoff A ~ 1 GeV, which also sets the energy-momentum scale for 
the validity of the model. That is, the model cannot be used for studying phenomena in- 
volving momenta running in loops larger than A. In particular, the model cannot be used 
to study quark matter at high densities pb ~ kp with kp > A, where kp is the quark Fermi 
momentum. One of the symptoms of this problem is the prediction of vanishing supercon- 



ducting gaps at high baryon densities, a feature of the model that is so 
use of a regulari zing momentum cutoff A of the divergent loop integrals 



ely caused by the 
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Recently, we |lj] have presented an alternative to the cutoff regularization of the NJL 
model, in which the one-loop integrals are reorganized through mathematical identities in 
a way finite integrals become separated from density-independent divergent integrals. The 
finite integrals are integrated without imposing any restriction on the integration momenta 
and the divergent integrals are related to physical quantities at the dynamical chiral sym- 
metry breaking in vacuum. As a result, instead of a vanishing gap like in the cutoff regular- 
ization one obtains a finite gap that grows with density. However, this does not mean that 
the model can be used at arbitrarily high densities. The point is that any implicit regular- 
ized integral contains an implicit scale, the regularization scale. When such an integral is 
fitted to a physical quantity in vacuum, we are implicitly also fitting this scale. Now, in a 
loop integral at nonzero density and/or temperature, the vacuum part is determined by the 
regularization scale and the finite integrals are governed by the chemical potential and tem- 
perature. When these scales are very different, one obtains nonsense results. Nevertheless, 
this mismatch of scales is less problematic in the implicit regularization scheme than in any 
other cutoff regularization scheme, in the sense that one can go to much higher densities 
and temperatures before the problem shows up. 



The root of the problem is of course the nonrenormahzabihty of the modeL In a renormal- 
izable theory all the dependence on the cutoff can be removed in favor of running physical 
parameters, like the coupling constants of QED and QCD. The running is given by the 
renormalization group equations of the theory and is controlled by an energy scale that is 
adjusted to the scale of the experimental conditions under consideration. In a recent pub- 
lication, Casalbuoni et al. |ll!] have introduced the concept of a running coupling constant 
for the NJL model to extend the applicability of the model to high density. The arguments 
are based on making the cutoff density dependent, using an analogy with an ordinary solid. 
In an ordinary solid, there is a natural maximum phonon frequency - the Debye frequency - 
that increases as the lattice spacing gets smaller. In quark matter, as the density of quarks 
increases, the quarks get closer together and therefore the ultraviolet cutoff A of the NJL 
model should be allowed to increase correspondingly. The cutoff can be changed consis- 
tently without spoiling the predictions of the model for the chiral properties of the vacuum 
if the four-fermion coupling constant G of the model is allowed to change with the cutoff. In 



Ref. [Uj], the cutoff dependence of G is obtained from the joint consideration of the divergent 
expressions for the pion decay constant f-,, and the gap equation for the constituent quark 
mass M : fixing the value of /tt at 93 MeV, the divergent expression of fj^ (in the chiral limit), 
regulated with a cutoff A, leads to a constituent quark mass M that is A dependent. The 
use of this M = M(A) into the gap equation gives rise to a coupling G that runs with A. 
Moreover, the running of the coupling is such that G{A) decreases with increasing A. This 
is certainly physically motivated and is also in accord with the interpretation [2, ll5| that 
the cutoff in the NJL model simulates - albeit in a crude way - the property of asymptotic 
freedom of QCD, in the sense that the coupling between quarks decreases as higher momen- 
tum scales are probed. A similar suggestion was made by Shakin and collaborators [l6| in 
studies of hadronic correlation functions. They found that results in qualitative accord with 
lattice QCD calculations can be obtained once a significant temperature dependence of the 
cutoff and of the coupling G is allowed. 

The aim of the present paper is to show that it is possible to extend the applicability 
of the NJL model to high densities and temperatures without the use of an explicit regu- 
larization of divergences. The basic motivation for avoiding an explicit regularization such 
as those commonly used like the three- or four-momentum cutoff, Pauli-Villars and proper- 
time regularizations, is that these lead to global and gauge symmetry violations, and to the 



breaking of causality and unitarity. Although in many situations these problems do not have 
great influences on the final numerical results, there are situations, however, where they do 
have drastic consequences, like in studies of correlation functions. Arguments have been 
used and tricks invented to circumvent such problems, but no satisfactory solution has been 
found - for a good discussion on these issues, see for example Refs. 



The implicit regularization scheme has been originally proposed in Ref. [201 and was used 



in different contexts 



21 



22|, including applications to the NJL model [ij, l23|, uS, l25| 



At the one-loop approximation, there appears only two divergent integrals, one is quadrat- 
ically divergent and the other is logarithmically divergent. Once the divergent integrals 
are assumed to be implicitly regulated, i.e. regulated without the specification an exphcit 
regulator, they will depend implicitly on a momentum scale that we denote by A. These 
integrals satisfy well-defined scaling relations, in that the integrals at some mass scale can 
be related to a combination of the same integrals at some another, arbitrary scale. These 
scaling relations allow to express the divergent parts of the amplitudes at finite temperature 
and density in terms of their counterparts at zero temperature and density. Moreover, fixing 
/,r in the vacuum, one can derive a scaling relation for the four-fermion coupling G which 
gives a running G{A) similar to that found in Ref. llj. All this is achieved, we reiterate, 
without specifying an explicit regularization scheme. 

In the next section we review the main aspects of the implicit regularization scheme \l4 . 
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251] that are most relevant for the present paper. We start specifying the 



Lagrangian density we are going to use and then explain the scaling relations and their use 
for obtaining the solution of the gap equation within this scheme. In order to make contact 
with results of the literature, we solve numerically this gap equation at finite temperature 
and density and show the results for the quark condensate as a function of temperature 
T and quark chemical potential fi. In Section llllt we make use of the scaling relations to 
obtain the running of the coupling G of the model and compare results with Ref. HI]. In 
Section IIVI we illustrate the use of our results for obtaining the critical temperature as a 
function of quark chemical potential for the spin-0 two-fiavor color superconducting (2SC) 
gap. Our conclusions and outlook are presented in Section |Vl 



II. CALCULATIONAL STRATEGY TO HANDLE THE DIVERGENCES 

For the purposes of the present work, it is sufficient to consider the simplest version of 
the model specified by the two-flavor SU{2) Lagrangian density 

Cnjl = i){%9- mo)^ + G{{i)i)f - (^7'^V')']- (1) 

Here ■?/' is the quark field operator (with color and flavor indices suppressed), mo is the 
current-quark mass matrix. At the one-loop approximation, the gap equation at finite 
temperature T and quark chemical potential /i, is given by (for simplicity we work in the 
chiral limit vtiq = 0) 

M = A8GMTy f-^- ^ — — -, (2) 

^J (27r)3 (z^„ + /i)2 + E(A;)2' ^ ^ 



with E{k) = ^/WTJP, and u;„ = (2n + 1)ttT, n = 0, ±1,±2,--- being the Matsubara 
frequencies. Performing the sum, one obtains 

M = A8GM [tI,uad{M) - /(T, /x)] , (3) 

where 



(2^4 ^2 _ ]Vf 2 

and 



quad[M ) — j ^^_^^ ,^^ ,_^^, (4j 



^^^'^^ = 7 (2^0^ 2E{k) ' ^^ 

with n±{k) being the quark and antiquark Fermi-Dirac distributions 

While I{T,^) is finite, IquadiM"^) is divergent. The divergent integral is an element of the 



systematization adopted in the calculational strategy [26| which we use in the present work. 
In 4D one-loop calculations, within the context of theories and models where the divergence 
degree is not higher than cubic, which is the case of the NJL model, besides Iquad{M'^) there 
also appears the basic logarithmically divergent object 

Iiog{M') = J (2,r)4(F-M2)2- (^) 

In the NJL model, IquadiM"^) is related to the constituent quark mass, as shown in Eq. ([3]), 
while Iiog{M^) is related to the pion decay constant in vacuum - see Eq. (I22l) . 



The traditional procedure to treat the divergent aniphtudes in nonrenormahzable models 
is the introduction of a regularization distribution in all loop integrals in order to render 
them convergent. In the context of the NJL model, with rare exceptions - such as in 



Refs. 



28 



29 



301 ] - the finite integral I{T, fis) containing the Fermi-Dirac distributions is 
also regularized. Due to the nonrenormahzable character of the model, the regularization in- 
troduced in this way cannot be removed along the calculations and therefore the predictions 
become associated with the chosen regularization. In practice, the amplitudes become func- 
tions of the regularization distribution parameters and consequently the energy-momentum 
dependence of the amplitudes, as emerging from the Feynman rules, are modified. As it is 
well-known, this leads to violation of fundamental symmetries like those associated to space- 
time homogeneity or to the introduction of non-physical thresholds that break unitarity. 
Having this in mind, an alternative strategy to handle the divergences typical of the 



perturbative calculations in QFT has been developed in Ref. [20|. The main idea is to 
avoid the critical step of the explicit evaluation of a divergent integral. The finite and 
divergent parts are separated by using a convenient representation for the propagators, in 
that all the dependence on the physical momenta is located in finite integrals which are not 
affected by regularization, thus avoiding the contamination with non-physical behaviors. 
The divergent parts obtained in this way are reduced to only a small number of basic 
objects after the adoption of a set of properties for purely divergent integrals, denominated 
consistency relations (CR), which are dictated by the requirements of symmetry preservation 
and elimination of ambiguities [7|]. The remaining divergent objects need not be solved 
explicitly; in renormalizable theories they are completely absorbed in renormalization of 
physical parameters, and in nonrenormahzable models they can be directly adjusted to 
physical quantities. 

There are many advantages in the adoption of the procedure just described, but two 
of them can be immediately noted: the treatment for the divergences is universal and the 
energy-momentum dependence of the amplitudes is not modified. Besides, a connection 
with traditional procedures can be made. In particular, all results can be mapped to the 
ones obtained within the context of dimensional regularization (DR), in situations where 
this method applies. On the other hand, the connection with regularizations based on reg- 
ularizing distributions requires some care in the case of nonrenormahzable models, due to 
the fact that the effects of the regularization on finite integrals are removed in our strat- 



egy. In this sense it is important to say that the prescription described above, which has 



been apphed in many problems involving divergences in QFT 2l|, represents a particular 
prescription to handle the amplitudes having divergences. Therefore, it should be consid- 
ered a procedure to make predictions with nonrenormalizable models in the same way as 
with any other regularization prescription. However, the important point to be noted is 
that, by construction many of the problems intrinsic to the traditional methods have been 
circumvented, like violations of symmetry relations among Green's functions, ambiguities, 
non-physical thresholds and scale invariance breaking. The convenience or usefulness of the 
method for the treatment of nonrenormalizable models should be judged by the consistency 
of its phenomenological implications. 

At one loop order, only two divergent integrals, Iguad{M^) and Iiog{M^) need be elimi- 
nated by some phenomenological adjustment. We first note that they are not independent 



quantities, they obey two scaling relations [27 1 

^quad 



{M') = /,w(A') + (m' - A^ Ii,,{\') 



+ 



(477) 



M^-A^-M^logj"^' 

2x i - ^M^^ 



2 



(8) 



Iio,{M') = Iio,{\') - (^log y-^j , (9) 

where A^ is the arbitrary mass scale adopted in the separation of finite and divergent parts. 
These relations allow us to relate the basic divergent quantities in two different mass scales, 
as we shall see in a moment. The imposition of scale independence implies in two properties 
for these objects 

dlquadi^ ) 



2 



dX 

dilog (A^) I 1 



ho, (A^) , (10) 



9A2 (47r)2A2' ^^^^ 

It is important to note that only finite quantities have been differentiated. The above prop- 
erties work like requirements that must be imposed in regularizations in order to maintain 
the scaling properties of the amplitudes. The above referred properties can be understood 
within the context of regularizations. For this purpose, we assume that each integral is 
regularized through an unspecified distribution f{k/A), where A is a parameter with the 
dimensions of momentum such that 

2, f d^k f{k/A) 



IquadiM ) - / (2^)4 p_M2' ^^^^ 



., , d^k f{k/A) 



Differentiating IquadiM"^), relation (ITOl) follows immediately. By differentiating (IT3l) . on the 

other hand, we obtain 

dliog (M2) _rd^k fik/A) 



dM^ "J (27r)4(A;2-M2)3- (^^^ 

Since the integral is finite and assuming the existence of the limit a^oo/(^/^) — ^i o^^ can 
extract the regulating distribution and perform the integration to obtain fITT]) . This last 
point illustrates the basic difference between our prescription and the traditional ones. 

The fact that the parameter A plays the role of a cutoff can be made clear even in the 
case of an implicit regularization. That is, it is possible to make the A dependence explicit 
without specifying the regulating function. This can be done by defining the dimensionless 
ratio Ma = M/A and writing Iquad{M'^) and Iiog{M'^) in terms of dimensionless integrals 
Jquad{Ml) and Jiog{Ml) as 

huad{M^) = A'jquadiMl), (15) 

hoqiM') = JioAMl), (16) 

with 

Jquad{Ml) = / -—-/l^, (17) 



(27r)4M2_M 



Jio,{Ml) = / ,,^,,'\[,,, , (18) 



. (27r)4(M2-Ml)2^ 

where the integration variable u is dimensionless. The integrals Jquad{Ml) and Jiog{Ml) 
obey scaling relations completely similar to that (IE]) and ([9]), as they should. A general 
parametrization for these integrals, obeying the properties ( ITOl) and ( ITTl) . can be constructed 
and written as 

t {Anf JiogiMl) = In (m1) + c,, (19) 

t (47r)' JquadiMl) = Ml In (m1) - M^ + c^. (20) 

It should be clear that if a regularization distribution does not furnish the above general form, 
the scaling properties of the amplitudes can be broken. In principle, different (consistent) 
regularizations should differ only by the values of the constants ci and C2- 



After this review on the general procedure to handle the divergences, let us turn our 
attention to its application to the NJL model. First, we have that the vacuum quark 
condensate is related to the scalar one-point function, which can be written as 

(V^V^)o = -12Mo iIquad{Ml), (21) 

where Mq is the mass obtained by solving the gap equation in vacuum. On the other hand, 
the pion decay constant f-,, can be related to the axial-pseudoscalar two-point function, such 
that 

fl = -l2Mltk,g{Ml). (22) 

The point now is that one can express the temperature- and density- dependent gap equa- 
tion ([3]) for quark mass M in terms of these vacuum quantities. This is so because in general 
the scaling relations ([8]) and ([9]) can be used to isolate vacuum contributions from tempera- 
ture - and density-dependent contributions in expressions for physical quantities, that is, it 
is possible to write IguadiM"^) and Iiogi^M"^) in terms of Iquad{MQ) and Iioq^Mq). Using the 
scaling relation in ((31), we get 



{Any 



I{T, /i) 



\ 12Mo ^ °^12Mo2 

1 



M^ - Mo^-M^ log (^ 



(23) 



where we have made use of Eqs. (1211) and (1221) . The corresponding expression for the case 
rriQ 7^ is a little more complicated than the one shown in Eq. fl23|) . since the expression for 
/^ contains, in addition to Iiog, a finite integral involving the pion mass. 

For illustrative purposes, we present in Fig. [1] the quark condensate as a function of 
temperature T and quark chemical potential fi - the results here are obtained solving the 
gap equation with rriQ = 5.5 MeV, G = 4.9 x 10"^ MeV~^ and the fitting mass Mq = 312 
MeV. As seen in Fig. [H this implicit regularization scheme gives the expected reduction of 
the quark condensate in medium. In addition, this reduction is in qualitative agreement 
with the result obtained with traditional cutoff regularization schemes - see for example 



Q 



Fig. 26 of Ref. [^. In the chiral limit, the order of the phase transition is first order. It is 



also worth mentioning that wit 
for the 7T- and a- meson masses 25 1 



1 these same parameters one is able to obtain good values 
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FIG. 1: One flavor quark condensate (ipip) as function of temperature T and quark chemical 
potential /i. 

Of course, the point of the exercise here was not only to show that this different way of 
handling the divergent integrals gives results in accord with the expected phenomenology. 
What should be noted here is that once we have eliminated Iquad^M"^) and Iiogi^M"^) in favor 
of {ipip)o and /^ we have kept A implicitly fixed. If one wants to change A without changing 
the low-energy results G has to run with A, as in Ref. Uj. The point here is that one can 
implement this using the scaling relations ([H]) and ([9]), which are independent of an explicit 
regularization. This will be done in the next Section. 



III. RUNNING OF THE COUPLING 



Fixing /tt = 93 MeV and using the expression of f-„ in terms of Iiog at the scales (M, A) 
and (Mq, Aq), one obtains M = M(A) as the solution of the transcendental equation 



M' 



,.2 A' 
Mq -^ exp 



3 M2 \M^ 



(24) 



Using the value of (V'V')o at the reference scale (Mo,Ao), and using the scaling relation 
to obtain Iquad at (M, A) in the gap equation, one obtains G = G{A) as 

-1 -(^^)o fM' M^\ fl 



48A"G(A) 



l2MaKl 



o^^o 



(47r)^ 



M2 
A^ 



A2 
Ml 



K 
A^ 



12MI 



A2 

^^0 



M2 / M^Ag ^ 
A^ ^^l^A2Mf^ 



(25) 
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where it is understood that M = M(A), as given by Eq. (12^ . 

This is our main resuh in the present paper. It is an interesting resuh, in that it was 
obtained without specifying any exphcit regularization, only very general scaling relations 
of the divergent integrals were used. In this sense, the result seems very general and robust. 
Once a temperature and density dependence for A is specified, G becomes also temperature 
and density dependent. 



350 



300 

> 

S 250 



200 



150 



1000 



1100 1200 

A (MeV) 



1300 



1400 



FIG. 2: Constituent quark mass M as a function of A. 



In Fig. [2] we present the numerical results for the A dependence of M, as obtained from 
the solution of Eq. (124|) after solving the gap equation. The constituent mass decreases as A 
is increased. If G were kept fixed, M would increase of course. But G decreases with A, as 
shown in Fig. [3], and the net effect is that M decreases. These results for M = M(A) and 
G = G{A) are in qualitative agreement with the results of Ref. 11|. 



The specification of a density and temperature dependence for A might seem a little 
arbitrary. But this needs not be the case, because the entire line of arguments can be turned 
around, instead of fixing /^ at some value, one could postulate in a physically motivated way 
a running behavior for G = G{A) and work backwards. What would change in this case? 
First, using the postulated G{A) in the gap equation, one would obtain a corresponding 
M(A). When this M(A) is replaced in the expression for f^^, in general /^r will be also 
A dependent, but this A dependence would be very weak since the integral for /^ is only 
logarithmically divergent. In this way, the chiral physics in vacuum would be maintained. 
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FIG. 3: Running of the coupling G(A). 

This is very interesting, since one could use the predicted density and temperature running of 
the QCD couphng constant for G and in this way fit the density and temperature dependence 
of A. In this way, one would be modeling, admittedly in a crude way, the asymptotic freedom 
of QCD in the NJL model. 

In closing this Section we reiterate that the purpose for making G to run with A is to 
extend the applicability of the model to high densities and temperatures. At high densities 
and temperatures, high momentum components are present in the system and a cutoff of 
the order of the chiral symmetry breaking scale invalidates the use of the model in such 
situations. In order to illustrate the use of the extension in practice we consider in the next 
Section color superconductivity in high-density quark matter. 

IV. COLOR SUPERCONDUCTIVITY 



For reviews and a cornprehensive list o 



ity see Refs. [31 
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34 



35 



36 



37 



' references on the subject of color super conduct iv- 
. Since the aim here is to illustrate the formalism 
we simplify matters by using the same Lagrangian density as used above, although for ob- 
taining a better phenomenological description of both chiral symmetry breaking and color 



13, MM. At 



superconductivity a more general four-fermion Lagrangian should be used 

the one-loop level the self-consistent equation for the superconducting gap is given by (we 

use the same letter G to denote both the diquark coupling here and the quark-antiquark 



13 



coupling in the last section) 
1 = 16^ 



d^k 



{2^lYkl-{k + ^ly-^^ 



+ (/i ^ -/i). 



(26) 



The integrals above are divergent. The application of the implicit regularization scheme to 
this problem proceeds as follows [ij]. Instead of introducing a cutoff in the integrals, the 
integrands are assumed to be implicit regularized and then manipulated in a way divergences 
are isolated in /i- independent divergent integrals. These divergent integrals can be related to 
the divergent integrals Iquad and Iiog of the problem of chiral symmetry breaking in vacuum 
though the use of the scaling relations discussed above. From the manipulation of the 
integrand in Eq. fl26l) also result finite integrals, and these are integrated without imposing 
any restriction to their integrands. Initially we consider T = 0. In this case the equation for 
the superconducting gap is given by 



1 



2G'(A)A^ 

A2 



4 {i,i,) 







3 M.kl 



TT^ 



Ml 

M 



A2 



4 
3 



'A2 
A^ 



Ml 



Jl 
Mi 



log 



[A^Al 



A2 



2 



A2 ^ \A^Mi 
/A^An 



'A2 



3Mi 7r2 ^\A^Mi) 



(27) 



where we used the manipulations of the integrand in Eq. (126|1 as explained above and shown 
explicitly in Ref. 14l |. 

One very interesting result of the application of the implicit regularization to the problem 
of color superconductivity is that, as shown with greater detail in Ref. p^, the supercon- 
ducting gap as a function of fi does not vanish at /i ^ A, as happens with the traditional 
cutoff schemes. The numerical result for A is shown by the solid line in Fig. |H This is 
interesting because a nonvanishing gap at high quark densities is predicted by QCD 41 1 



see also the reviews in Ref. 
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35 



36 



371 . |38| . Of course, if the implicit 



regularization scale A is kept fixed, the use of the model at high densities is questionable. 
However, we are able to make explicit the A-dependence in the gap equation (and in other 
physical quantities as well) by extracting the A dependence from the implicit regularization 
function. In this way, one can very easily extend the applicability of the model to larger 
values of /i by allowing a running G(A) and a /i dependence for A. 
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FIG. 4: Zero temperature superconducting gap as a function of quark chemical potential for three 
different values of a ii □ 
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FIG. 5: Critical temperature for the superconducting gap as a function of /x for different values of 
a in Eq. (l28]l . 



For the density dependence of A, as said above, there is a great deal of arbitrariness. As 
discussed, one could fix this dependence, for instance, by matching the density dependence 
of G'(A)A^ with the prediction of perturbative QCD for the running of the QCD coupling 
constant a^ at high densities - the one loop prediction is that the coupling decreases log- 
arithmically with fi for large values of fi. However, for our purposes here of showing the 
qualitative results only, we use the simple formula for /i > /io = 235 MeV, 



A = An 



1 + a log 






(28) 
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where a is a constant. When this is used into the expression for G'(A)A^, one obtains that 
G(A)A^ decreases logarithmically with fi, mocking up in a rather crude way the prediction 
of perturbative QCD for the running of a^ with fi. In Fig. |l]we show the numerical results 
A as a function /i for different values of a. We use for the diquark pairing strength the value 
G = 3.1 GeV-^. The other parameters are Mq = 312 MeV and Aq = 932 MeV, which are 
the values obtained in vacuum for /^r = 93 MeV and {ipil')o = (~ 250. MeV)^. The results 
show the expected behavior of the gap growing faster with yU as A increases with fi, i.e. as 
a increases. 

For completeness, we include the effects of temperature. We calculate the critical tem- 
perature Tc above which A = 0, for different values of fi. We rewrite Eq. (1261) for A = 
including the effect of temperature 



1 = 16G 



T — o tanh \ —\k 



fi\ 



(29) 



where 



tanh 



'P 



k±fi\) = l-2n±(/c) 



(30) 



with 



< (k) 



(31) 



Using the same manipulations to isolate the divergent integrals as explained before and after 
some algebraic effort we can write 



— = A^< 
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3 I Ag ^ A2 



J IT 



3 MnA2 
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(33) 
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with 

In Fig. ([5]) we plot the critical temperature as function of /i. In the solid line we have 
kept the coupling G fixed and in the dashed and dotted lines the coupling is a function of /i 
through the /i dependence of A as in Eq. ( l28l) . As seen, the critical temperature increases 
with /i, and this increase is faster as a increases for large /i. The values of T^ obtained within 
this scheme seem a little larger than obtained with cutoff regularization. 

V. CONCLUSIONS AND PERSPECTIVES 

We have considered an extension of the NJL model to high densities and temperatures as 
proposed in Ref. [11|] and used an alternative strategy to handle ultraviolet divergences. This 
extension is implemented by allowing the regularization scale A to increase at high densities 
with the simultaneous decrease of the coupling G. Making use of the scaling relations of 
Eqs. (IHD and (j9]), and the definitions in Eqs. ( TT5l) and ( TT6l) . the two one-loop divergent 
integrals Iquad and Iiog at scales Aq and Mq can be related to the same Iquad and Iiog at 
different scales A and M. 

Although our numerical results are in qualitative agreement with the ones of Casalbuoni 
et al. [11] for the yU dependence of the superconducting gap A, we believe the present 
approach is different from their results in several aspects. Perhaps the most important 
one is the fact that no specific regularization distribution was used to obtain the running 
of G with A. Proceeding this way it is possible to circumvent many troubles commonly 
found when traditional regularization schemes are used, like those based on a three- or four- 
momentum cutoff, or Pauli-Villars and proper-time regularizations, lead in general to global 
and gauge symmetry violations, and to the breaking of causality and unitarity. Although 
in some circumstances such aspects might not be of crucial importance, in many others. 



like in hadronic correlation functions, remarkable differences appear [42|. Among them, 
only physical thresholds (independent of A) are present allowing thus the preservation of 
unitarity. In addition, it is easy to show that causality is also preserved by checking that 
the amplitudes obey the correct dispersion relations. Gauge symmetry is also preserved and 



the ambiguities are eliminated at the one-loop level, as shown in Ref. 23| in the context of 



the gauged NJL model. A complete discussion of these matters are discussed in detail in a 

17 



separate publication |25 |. 



Another interesting aspect of the present approach is that it can be used with heavy 
quarks. The scahng relations involve the ratio of the quark mass to the implicit regulation 
scale and naturally take into account short-distance effects as the quark mass increases. This 
is very important in connection with studies of heavy-quark bound states in highly excited 
quark matter. 

Acknowledgments: Work partially supported by the Brazilian agencies CAPES, CNPq, 
FAFERJ and FAPESP. 
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